NEW MEASUREMENT OF METAL ABUNDANCE IN THE 
ELLIPTICAL GALAXY NGC 4636 WITH ASCA 
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ABSTRACT 



High quality X-ray spectra of NGC 4636 are obtained with ASCA. 
Theoretical models are found unable to reproduce the data in the Fe-L line 
region. Spectral data above 1.4 keV indicate that Mg to Si abundance ratio 
is ~ 1 solar. Assuming that the abundance ratios among a-elements are the 
same with the solar ratios, spectral fit with increased systematic error in the 
0.4-1.6 keV range gives abundances of ct-elements and Fe to be both ~ 1 solar 
by taking one solar of Fe to be 3.24x10"^ by number. These new abundance 
results solve discrepancy between stellar and hot-gas mctallicity, but still a low 
supernova rate is implied. We also detect strong abundance gradients for both 
a-elements and iron in a similar fashion. The abundance is ~ 1 solar within 4', 
and decreases outerwards down to 0.2 ~ 0.3 solar at 10' from the galaxy center. 
Dilution due to an extended hot gas is suggested. 



Subject headings: galaxy: abundance — X-rays:galaxies — galaxy individual 
(NGC 4636) 
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1. Introduction 

The interstellar medium (ISM) of elliptical galaxies is thought to be considered to be 
an accumulation of stellar mass loss and supernova ejecta. While standard supernova rates 



predict the ISM to have a high metallicity of several times the solar value (e.g. [Loewenstein 



|et al. 1991| ; [Renzini et al. 1993| ), the observed ISM metallicity in fact fell below half a 



solar, even lower than the stellar metallicity (e.g. |Awaki et al. 199^ ; [Mushotzky et al" 



[Matsushita et al. 1994| ; [Arimoto et al. 19971 ; |Matsumoto et al. 19970 . This strong 



discrepancy calls into question our current understanding of supernova enrichment and 
chemical evolution of galaxies. 

NGC 4636 is one of the relatively isolated yet most luminous elliptical galaxies, 
both in optical (Lb=2.9 x IO^^Lq assuming a distance of 17 Mpc) and in X-ray {Lx= 
3.8xl0^^erg s~^; Einstein ) band (e.g. [Forman et al. 1985| ; [Fabbiano et al. 1992| ). Using 



ROSAT, Trinchieri et al. (1994) discovered a very extended X-ray emission surrounding 
this galaxy, out to ~ 18'. The first ASCA observation of NGC 4636 in the PV (Performance 
Verification) phase yielded the very low ISM metallicity ( [Awaki et al. 199^) , together with 
abundance and temperature gradients ( pVIushotzky et al. 1994|) . 



We have re-observed NGC 4636 with ASCA ( [Tanaka et al. 19"9^) for an extremely long 
time; over 200 ks. This has allowed Matsushita (1997) to perform a much deeper study of 
the extended X-ray component after ROSAT. In this paper, we utilize the overall ASCA 
data including this long exposure, to study the spectral properties of NGC 4636 and look 
into the abundance problem. 
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2. Observations 

NGC 4636 has been observed from ASCA twice. Following the PV observation (1993 
July 22, with the SIS in 4CCD mode), the second much longer observation was conducted 
in the AO-4 phase from 1995 December 28 through 1996 January 4, with the SIS in the 
ICCD mode. We discarded the data taken under cut-off rigidities less than 6 GeV c~^, or 
elevation angle less than 5° and 20° from night and day earth respectively. This has yielded 
exposure times of 36 ks (with the GIS) and 39 ks (with the SIS) for the PV observation, 
and those for the AO-4 observation are 172 ks (GIS) and 215 ks (SIS). 



We accumulated on-source spectra within 4 times the effective radius, rg = 1'.7 ( Paber 



et al. 19"5PD, centered on NGC 4636. Since the SIS response has changed with time, we 



treat the SIS spectra from the two observations separately. The background spectrum was 
obtained by integrating the blank-sky data over the same region of the detector. Figure |l| 
shows the background-subtracted SIS spectrum for the AO-4 data. 



3. Results 

We jointly fit the two SIS (PV and AO-4) spectra and one GIS spectrum, with a 
standard two component model ( [Awaki et al. 19^ ; Matsushita et al. 1994 ; Matsumoto 



et al. 1997| ). The model consists of a thin thermal emission from the ISM with free 



temperature kT and free metallicity, and a thermal bremsstrahlung with temperature fixed 
at 10 keV representing the contribution from low-mass X-ray binaries. Both components 
are subjected to a common interstellar absorption A^h- Although the GIS ( |Ohashi et al. 



19961 ; [Makishima et al. 1996| ) is less sensitive than the SIS to the low-energy (e.g. < 1 keV) 



atomic lines, it can constrain the hard component better than the SIS, and is fully sensitive 
to the Si-K and S-K lines. Therefore, the joint use of the two instruments is essential. In 
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this paper we adopt for the solar iron abundance the 'meteoritic' value, Fe/H= 3.24 x 10 
by number ( [Anders fc Grevesse 1989| ) 



3.1. Problems with the spectral analysis 

As the first-cut spectral analysis, we represented the ISM component by the plasma 
emission model of Raymond & Smith (1977; hereafter R-S model) with solar abundance 
ratios. The best-fit model parameters turned out to be consistent with those of Awaki et 
al. (1994), but the fit was totally unacceptable (Table 1). We then allowed the abundances 
to deviate from the solar ratios, by dividing heavy elements into two groups to estimate 
relative contributions from type la and type II supernovae (SNe). One group consists of 
so-called a-elements, O, Ne, Mg, Si and S, which are assumed to have a common abundance 
Aa- The other group includes Fe and Ni, with a common abundance Ap^- Abundance of He 
is fixed to be 1 solar. The abundances of the other elements are assumed to be the same as 
a-elements, although their effect on 1 keV spectrum is negligible. The fit incorporating 6 
parameters {Ac, Ap^, kT, N^, and two normalizations) still remained far from acceptable, 
as shown in Table 1. We further replaced the R-S model with MEKA model ( |Mewe et al.| 
T985t IMewe et al. 1986| ; |Kaastra 199'^) , MEKAL model ( |Liedahl et al. 1995|) , or Masai 



model ( [Masai 1984| ), but none of them were successful. Multi-temperature ISM models did 



not improve the fit, either. 

In these fits, the data-to-model discrepancy is always largest around the Fe-L complex 
region (0.8 ~ 1.4 keV). In addition, the fit residuals and the derived physical quantities both 
depend significantly on the plasma emission codes. These reconfirm the serious problems in 
the theoretical Fe-L line modeling, as pointed out by Fabian et al. (1994) and Arimoto et 
al. (1997). Furthermore, we have found strong false couplings between Ac and Ape, arising 
from the following two effects. On one hand, the fitting algorithm tries to reduce the Fe-L 
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discrepancy, by adjusting intensities of the 0-K and Ne-K lines which overlap the Fe-L 
complex; this strongly affects A^, since the data have the highest statistics in this energy 
range. In turn affects A-pe, because the bound-free emission from oxygen and neon acts 
nearly as a continuum to the Fe-L lines and controls their equivalent widths. These effects 
make both Ap^ and A^ highly unreliable. 



3.2. Abundance of a-elements 

In order to avoid these problems, we tentatively restricted the energy range for the SIS 
spectral fit to > 1.4 keV which is unaffected by the Fe-L, 0-K and Ne-K lines. For the GIS 
spectrum, we used an energy range above 1.6 keV, because of its poorer energy resolution. 
We re-arranged heavy elements into 4 groups; Mg; Si; S; and other elements including 
O, Ne, Fe and Ni (hereafter Aothors)- Elements in the last group contribute to the > 1.4 
keV spectrum mainly via free-bound continua (Fe-K and Ni-K lines are undetectable). 
Then, the fit has become acceptable with the reduced of 1-06 for 220 dof.. The derived 
temperature, 0.73 keV (0.67-0.76 keV with 90% confidence for a single parameter), agrees 
with those derived using the whole energy range. We have further found that the Mg to Si 
abundance ratio always agrees with the solar ratio (Fig. 2), even though their abundances 
correlate considerably with ^others- The S abundance may be similar, or slightly higher 
by about 30 %. On the other hand, the stellar and supernova nucleosysnthesis predicts 
O, Ne and Mg to follow the solar abundance ratios ( [Thielemann et al. 1996 ). Therefore, 



it is reasonable to assume that all the a-elements (O, Ne, Mg, Si and S) have the same 
abundance in solar units. 

We can now safely return to the original element grouping into A^ and Ap^,, and refit 
the spectra above 1.6 keV (SIS) and 1.7 keV (GIS). The fit is acceptable with the reduced 
of 1.06 for 210 dof., and the derived temperature is 0.75 (0.72-0.78) keV. Now, Aa is 
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determined relatively well (Fig. 3), because it relies upon the Si-K and S-K lines instead of 
using the unreliable 0-K and Ne-K lines that are now ignored. Unlike the full-range fit, 
Aa depends on Ape only weakly through Fe contribution to the free-bound continuum. 
These results unambiguously show that Aa is roughly one solar, unless Ape is extremely 
high. 



3.3. Iron abundance 

Our next task is to better constrain Ape- This requires us to utilize the Fe-L line 
information, but simply reviving the energy range below 1.6 keV would bring us back to the 
original problems described in §3.1. In order to relax the too strong constraints imposed 
by the data in the Fe-L region, we have decided to assign 20% systematic errors to each 
spectral data in the 0.4-1.6 keV range. This is because the data-to-model discrepancy and 



the model-to-model differences are both ~ 20 % around the Fe-L region ([Masai 1997|) . We 
have thus fitted the SIS/GIS spectra jointly with the two-component model, incorporating 
the increased systematic error and again employing the element grouping into A^ and Ape- 

The fit is acceptable with either of the plasma models (Table 1), and the confidence 
regions on the A^ vs. Ap^, plane for different plasma emission codes have come to a similar 
region (Fig.3). As a consequence, we can now conclude with confidence that A^ and 
Ape are both consistent with ~ 1 solar, within a factor of 2. In addition, the ratio between 
Aa to Ape coincides with the solar ratio within ~ 20%. 



3.4. Metal distribution in the ISM 



In order to study the metal distribution, we accumulated the GIS and SIS spectra 
in three concentric annular regions with radius range, 0-4', 4'-8', and 8'-12' centered on 
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the galaxy. We assume that the ISM is spherically symmetric. This "ring-sorted analysis" 
is however limited by the spectral mixing effect among different sky regions, due to the 
extended point- spread-function of the ASCA X-ray Telescope (XRT; [Serlemitsos et al- 
ii 995|) . This tends to smear out any radial gradient in the ISM properties. Following the 
analysis in §3.3, we include 20% systematic error in the 0.4-1.6 keV data. We fit only with 
the R-S model, since the difference due to plasma emission code is only several tens of % 
after inclusion of the systematic error. An acceptable fit is obtained in all three regions 



As shown in Table 2 and Fig.4, the derived Aa and Ape are ~ 1 solar within r < 4', 
and drop to 0.2 ~ 0.3 solar at r ~ 10'. This confirms the existence of strong abundance 
gradient in the ISM in NGC 4636, and reveals that both a-elements and Fe equally show 
the abundance gradient. Their ratio is close to the solar value at all radii. The absolute 
values of the abundances are higher than those derived in Mushotzky et al. (1994), who 
assumed a solar ratio in the metals, but the temperatures and hydrogen column densities 
are consistent. 



By carefully examining the abundance ratios and uncertainties in the Fe-L complex, 
we have concluded that the ISM abundances of Fe and a-elements in the X-ray luminous 
galaxy, NGC 4636, are in fact as high as about 1 solar within r < 4'. The abundance ratio 
between Aa and Apc is very close to the solar ratio within ~ 20%. This means that both 
Type la and Type II SNe ejecta make significant contribution in the ISM. 

The ISM iron abundance is expected to be given by 



(Table 2). 



4. 



Discussion 




(1) 
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( |Loewenstein et al. 1991 ; Renzini et al. 1993 ), where < > is the average iron abundance 



of the stars in units of the solar abundance, h^Q = Hq/50 the Hubble constant in units of 
50 km s~^ Mpc~^, "i^sN i is the present rate of SN la in ellipticals in units of the rate as 
estimated by Tammann (1982), and M^'^ is the iron yield per SN la event. The stellar Fe 
abundance, < Zj"^ >, averaged over the whole galaxy is 0.74 solar ( [Arimoto et al. 1997|) . 
Since our revised value of Ape is fairly close to these stellar values, the previously- quoted 
severe problem that the ISM abundance was lower than the stellar metallicity has been 
mostly removed. 

We can solve for the relative contribution from SN la and SN II on a-elements and Fe. 
If we assume that SN II products have 3 times higher value of Aa than Ap^ ( [i'hielemann 



et al. 1996| ) and SN la produces only Fe, the SN II contribution necessary to produce 



the observed A^ (1 solar) can enrich Fe to Ap^ = 1/3 solar. The rest of Ap^, i.e. 2/3 
of the best-fit Fe abundance, has to be explained by SN la. Using equation (1), we can 
estimate the SN la rate ^sm- In order to derive an upper limit for ^sm, let us assume that 
stars contain no Fe from SN la. This gives an equation only for the SN la contribution, 
■^KM.SNia > St^sN I (m^'^/O.TMq^ /igg. The left hand side is now estimated to be ~ 0.7 
solar, so that the condition for {}sm is "(^sni^so < 0.14. Therefore, the present abundance 
result suggests that SN la rate would be less than 1/7 of Tammann's rate for h^Q = 1. 
Considering the uncertainties in A^ and Ap^, this is consistent with the recent estimation 
of "i^sNi ~ 1/4 by Capperallo et al. (1993) based on the supernova counting. A systematic 
study of the ISM abundance involving many galaxies would be neccessary to further look 
into this problem. Also, Fe-L diagnostics need to be closely examined in order to reach 
consistent answers among the models ( [Arimoto et al. 1997 ). 



We have confirmed and refined the strong abundance gradient feature in the ISM 
discovered by Mushotzky et al. (1994). In addition, we have for the first time discovered 



strong abundance gradient in a-elements, which is very similar to that of Ap^- There 
may be an internal metallicity gradient in the ISM within the galaxy, as suggested from 
the significant gradient in the stellar metallicity within re ( Pavies et al. 19931) . However, 
no information is available about the stellar metallicity in larger scales. There exists an 
extended X-ray emission around NGC 4636 with a radius > 100 kpc, as discovered by 
Trinchieri et al. (1994), and confirmed with ASCA, by Matsushita (1997). Since the hot 
gas responsible for the extended emission has a very low metallicity as seen in the data in 
r=8-12', the abundance gradient can be explained in terms of the growing contribution 
from the metal poor gas in the outer region. The similarity in the slope of of Aa and 
A-pf, gradient suggests common dilution due to external gas. 



The authors are grateful to K. Matsuzaki who have made possible the exceptionally 
long observation from ASCA with incredibly efficient command arrangements. 
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Fig. 1.— The SIS spectrum (AO-4 observation) of NGC 4636 fitted witfi a two- 
component model, consisting of a soft R-S component (dotted lines) and a hard 
bremsstrahlung component with kT=10 keV (dotted line). Abundance of heavy elements 
are grouped into a-elements and iron group. The bottom panel shows data-to-model ratios 
for R-S model (crosses) and MEKAL model (diamonds). 

Fig. 2. — Two-parameter 90% confidence contours of abundances of Mg (thin hne), Si (bold 
line), and S (dashed hues), shown against the ^others- These were obtained from the joint fit 
to the spectra above 1.4 keV (SIS) and 1.6 keV (GIS). 

Fig. 3. — 90 % confidence contour of the correlation between Apc vs. Aa obtained by 
fitting the spectra above 1.6 keV (SIS) and 1.7 keV (GIS), using the R-S model (dot-dash), 
and those obtained by R-S (bold), MEKAL (dash), MEKA (thin) and Masai (dot) models 
fitting the whole energy band, including 20 % systematic errors in the energy range 0.4-1.6 
keV. 

Fig. 4. — Projected radial profiles of (crosses) and Ape (diamonds). 
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Code 


kT 






A 


A T 

Nh 


2 / 




(keV) 


(solar) 


(solar) 


(lO^'^cm 




R-S 


0.76 




0.48 




5 


1112/373 


R-b 


0.75 




0.54 


0.50 


4 


1093/372 


MhvKA 


0.59 




1.66 


1.49 


9 


1235/372 


MbKAL 


0.66 




0.76 


0.71 


10 


727.5/372 


Masai 


0.71 




1.50 


1.68 


1 


786.6/372 


Including 


systematic error in Fe-L region (0.4-1.6 keV) of the SIS spectra 


R-S 


0.76 


+0.02 
-0.02 


0.95 


U. '0 -0.09 


11? 


223.0/330 


MEKA 


0.63 


+0.02 
-0.02 


1.10 t',it 


-1 -19 +0.37 
^■^^ -0.20 




337.4/330 


MEKAL 


0.72 


+0.01 
-0.02 


1 iq +0.16 

^■^'^ -0.19 


0.88 




247.4/330 


Masai 


0.72 


+0.02 
-0.02 


1 1 7 +0.25 
-0.20 


-r qq +0.30 
l.OO _o.23 


9 +2 
^ -2 


271.5/330 



Table 1: Summary of joint GIS and SIS spectral fits for the NGC 4636 data accumulated 
within r < 4re (6'. 8). Errors show 90% confidence hmits for a single parameter. 



"The definition of solar iron abundance: Fe/H=3.24xl0 ^ (by number) 



r 


kT 






Nh 




arcmin) 


(keV) 


(solar) 


(solar) 


(1020cm-2) 




0-4 


0.75 


1.19 1^:1 


0.93 IgiJ 


9 +3 
^ -2 


182.5/264 


4-8 


0.81 l^i 


0.53 t',it 


0.48 IH? 


lit 


202.3/259 


8-12 


n 77 +0.05 


0.37 f'oil 


25 

^■^'J -0.08 




149.0/194 



Table 1: The fitting results for the ring- cut SIS spectra by adding 20 % systematic errors in 
0.4-1.6 keV. The Raymond-Smith model has been used. 
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